Quantum magnets with weakly confined spinons: 
Multiple length scales and quantum impurities 
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In magnets with strong quantum fluctuations, paramagnetic ground states with or without con- 
finement of spinon excitations can be realized. Here we discuss the physics of the confined phase in 
the vicinity of a confinement -deconfinement transition: this regime, likely relevant to a multitude 
of frustrated spin systems, is characterized by multiple length scales. In addition to the magnetic 
correlation length, a confinement length can be defined, which can be probed, e.g., by local static 
measurements near non-magnetic impurities. We illustrate the ideas by explicit calculations for 
dimerized spin chains, but our qualitative results remain valid in higher dimensions as well. In 
particular, we study crossover from weak to strong confinement visible in the antiferromagnetic 
polarization cloud around a non-magnetic impurity. We also discuss the effective magnetic interac- 
tion between impurities, relevant for impurity-induced magnetic order, and consequences for nuclear 
magnetic resonance and neutron scattering experiments. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee 
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I. INTRODUCTION 

The last decades have seen tremendous interest in ex- 
otic phases of correlated quantum matter characterized 
by elementary excitations with fractional quantum num- 
bers. Some examples are Luttinger liquids p spin liquids, 2 
and fractional quantum Hall systems.— While it has long 
been known that fractionalization generically occurs in 
one space dimension (ID), the situation in higher dimen- 
sions D is quite different. It has become clear that frac- 
tionalized phases in D > 2 can be connected to decon- 
fined phases of certain lattice gauge theories^ 

Quantum antiferromagnets have played a central role 
in investigations of exotic phases^ Such systems can ex- 
hibit both antiferromagnetic and paramagnetic ground 
states. In the latter situation, which is generic to ID 
due to the Mermin- Wagner theorem, quantum fluctua- 
tions prevent long-range order of the spins. In D > 2, 
quantum paramagnetic states may obtained by the for- 
mation of singlet bonds of spins 1/2. For lattices with 
one spin 1/2 per unit cell, the singlet bonds can form ei- 
ther crystalline order accompanied by lattice symmetry 
breaking, leading to so-called valence bond solid (VBS) 
phases, or short-range ordered liquid states without any 
broken symmetry.? Remarkably, such spin liquids can be 
shown to be fractionalized, with elementary excitations 
being spin- 1/2 spinons. 

The search for microscopic realizations of spin liquids 
in D > 2 remains difficult, because the relevant spin 
models display geometric frustration: this reduces the ap- 
plicability of otherwise powerful Quantum Monte Carlo 
techniques due to the so-called sign problem. The first 
unambiguous demonstration of a stable phase with frac- 
tionalization in two-dimensional (2D) magnet was for the 
triangular-lattice quantum dimer model, which features 
both VBS and gapped spin liquid phases, the latter dis- 
playing topological order^ Models with exact gapped 
spin-liquid ground states have been constructed, like the 



toric code model- and the spin- 1/2 honeycomb model^ 
both due to Kitacv. Hciscnberg models with longer-range 
frustrated interactions are believed to be good candidates 
for topological spin liquids as well, but with the exception 
of the triangular-lattice ring-exchange models no conclu- 
sive numerical results are available. (Gapless spin liquids 
with power-law correlation functions do not play a role 
in the subsequent discussion.) 

A special form of fractionalization has been argued to 
be possible at a continuous quantum phase transition 
between two conventional phases with different broken 
symmetries*^ A candidate for this scenario of deconfined 
quantum criticality is the Neel-VBS transition of the 
square-lattice Heisenberg antiferromagnet supplemented 
by a ring-exchange interaction*^ An interesting predic- 
tion of the deconfined critical theory is the existence of 
two distinct small energy scales on the VBS side of the 
quantum critical point: primarily, the transition is be- 
tween a Neel state and a spin liquid, and the VBS is a 
secondary instability of the spin liquid. Consequently, 
dimerization and confinement set in at a much smaller 
energy scale compared to paramagnetic behavior. 

A common feature of (i) a confined paramagnet near 
a confinement-deconfinement transition and (ii) a VBS 
phase near a deconfined critical point is that, in addi- 
tion to the magnetic correlation length £ mag , there exists 
a much larger length scale associated with confinement, 
which wc dub confinement length £ con f . In the following, 
we shall understand £ C onf as the typical spatial extension 
of a two-spinon bound state. Importantly, both £ con f 
and £ mag can be large compared to microscopic lengths, 
opening the possibility to find universal phenomena char- 
acterized by two length scales, in contrast to the conven- 
tional situation with £ mag being the only large length. In- 
deed, a recent numerical study-i 2 - of the frustrated regime 
of the J1-J2-J3 Heisenberg antiferromagnet on a square 
lattice gave evidence for two distinct length scales, with 

£conf ^ {mag ■ 
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The purpose of this paper is a detailed discussion of 
physical properties of quantum paramagnets near decon- 
finement. We believe this regime to be relevant for a 
rather broad class of frustrated magnets: Even if the el- 
ementary excitations are ultimately confined, there will 
be signatures of deconfincment on shorter length scales, 
which can be probed experimentally. The presence of two 
length scales, £ con f and £ m ag, strongly modifies a number 
of phenomena, both in statics and dynamics. We present 
general arguments and support them by explicit calcu- 
lations for ID Heisenberg models. Although the mech- 
anisms of fractionalization are quite different in ID and 
higher dimensions, we argue that various consequences 
for observables are similar, and hence our results remain 
qualitatively valid in D > 2 as well. 



A. Observables 

In the simplest picture, weak confinement can be un- 
derstood in terms of a (weak) linear potential, V(x) — 
a\x\, describing the effective interaction between two 
spinons. The two-particle problem then has only bound- 
state solutions, i.e., all excitations are of two-spinon type. 
The lowest two-spinon triplet state takes the role of the 
conventional spin-1 (triplon) excitation, its size defining 
the confinement length £ con f, with £ con f — > oo for a — > 0. 

A static measurement of £ C onf is possible, e.g., by spa- 
tial pinning of one of the two partners of the bound state. 
This can be achieved by doping with non-magnetic im- 
purities: removing one spin- 1/2 from the system breaks 
a singlet and liberates a spinon which now interacts with 
the vacancy (holon) . Provided that the binding between 
spinon and holon is identical to that between two spinons, 
the size of the spinon-holon bound state gives a measure 
of £conf • As we show below, the size of the antiferromag- 
netic polarization cloud around the vacancy is given by 
£conf in the weakly confined regime (whereas it is deter- 
mined by £mag in the conventional case of strong confine- 
ment). This polarization cloud, which can be detected 
e.g. by local magnetic probes, also determines the effec- 
tive magnetic interaction between impurity moments and 
is hence relevant for impurity- induced magnetic order. 

Weak confinement shows up as well in the dynam- 
ical properties, as it implies the existence of multiple 
bound states below the two- or three-triplon continuum. 
Hence, inelastic neutron scattering experiments, measur- 
ing the dynamic spin susceptibility, will detect a series of 
sharp dispersing modes below the threshold of the multi- 
particle continuum. 

With increasing confinement potential, £ con f decreases. 
The criterion for the crossover to strong confinement may 
be defined as £ C onf ~ £mag- This is expected to coincide 
with the condition that the energy spacing between the 
lowest two-spinon bound states becomes comparable to 
the spin gap, and hence the higher bound states are no 
longer discernible in the dynamical susceptibility. 



B. Outline 

The remainder of the paper is organized as follows. In 
Sees. HlHVl we discuss the physics of weak confinement in 
one dimension: Sec. HT1 summarizes the well-known prop- 
erties of the frustrated dimerized spin chain, described by 
the J1-J2S model, and the methods applied in the fol- 
lowing sections. Numerical results for the undoped spin 
chain are discussed in Sec. IIIII while analytical and nu- 
merical calculations for the system doped with a single 
vacancy are presented in Sec. [iVl Sec. [V] describes ex- 
perimental implications of weak confinement for nuclear 
magnetic resonance (NMR) and neutron scattering mea- 
surements and for impurity-induced magnetic order. Fi- 
nally, in Sec. I VII we discuss the relevance of our findings 
for higher dimensions. A brief summary of our results 
closes the paper. 



II. FRUSTRATED DIMERIZED SPIN CHAINS 

In one space dimension, spinon excitations are quite 
generically deconfined. The standard spin-1/2 Heisen- 
berg chain has a gapless excitation spectrum, while frus- 
tration can induce spontaneous dimerization leading to 
gapped (but still deconfined) spinons. To induce confine- 
ment, one has to resort to models with explicitly broken 
translational symmetry. 



A. Model and phase diagram 

We consider the so-called J1-J2-8 model, describing a 
frustrated dimerized antiferromagnetic (AF) spin chain: 

2N 2N 

H ch = J 1 ^(l-(-I) l 5)S. i -S 4+1 + JaJ^Si'St+a- 

»=i »=i 

(1) 

Here represents a spin-I/2 at site i, J\ and J2 are the 
antiferromagnetic first- and second-neighbor exchange 
couplings, and 5 denotes the degree of dimerization, see 
Fig. [Ha). (We will employ J\ = 1 in the following.) 

The ground-state phase diagram of the model fl} in 
the J2-8 plane is well known ) 13 ' 14 Fig. [21 and will be 
summarized in the following. We start with 5 = 0. For 
small J2 < J%Ci the ground state is unique and the ex- 
citation spectrum is gapless. At the critical value J^c — 
0.24I1 ) 15 ' 16 a continuous phase transition of Kosterlitz- 
Thouless type into a gapped phase takes place. Here, the 
system breaks an Ising symmetry accompanied by bond 
order, i.e., spontaneous dimerization. At J2 = 0.5, the 
so-called Majumdar-Ghosh (MG) point, the two ground- 
state wavefunctions are given by a nearest-neighbor sin- 
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FIG. 1: (color online) Schematic representations: (a) Model 
Hamiltonian ((TJ. The thick solid (black) lines correspond to 
a strong bond [Ji(l + <$)], the dashed lines to a weak bond 
[Ji(l — 5)], and the thin solid (red) lines to the J2 coupling, 
(b) State |0). The ellipses represent singlet states between 
neighboring sites, (c) State \nm). The two spins at sites 
2n — 1 and 2m are separated by a singlet pattern. 



glet product states^ 



|0) = [!,2][3,4][5,6]...[2iV-l,2iV], 



|0) = [2iV,l][2,3][4,5]...[2iV-2,2iV-l] 



(2) 



with denoting a normalized singlet combination of 
the spins at sites i and j, 



^0 mih 



lh\ T)j] 



(3) 



The state |0) is illustrated in Fig.rjjb). 

In the gapped phase, elementary excitations are topo- 
logical defects (solitons or spinons) separating the states 
|0) and |0). Shastry and Sutherland 18 performed vari- 
ational calculations for the MG point using a reduced 
Hilbert space {\nm)}, where the states |nm) consist of 
two spins at sites 2n — 1 and 2m separated by nearest- 
neighbor singlets, Fig. [IJc). It was shown that the low- 
lying excitations are pairs of propagating spinons, i.e., 
deconfined spin- 1/2 (fractionalized) quasiparticles. The 
upper bound for the spinon dispersion, = (l/2)(5/4 + 
cos2fc), was later confirmed by numerical calculations. 19 

A finite dimerization 5 opens a gap for all J 2 . The 
perfectly dimerized states |0) and |0) remain eigenstates 
of H c h along the so-called "disorder line" 2J 2 +5=1, 
but their degeneracy is lifted by 5 > 0. The energy 
of the ground state |0) is E = (-3/8) (1 + 5)(2N). It 
was also shown that the static structure factor S(q) = 
(1/2N) exp(i<7i?j)(So • Sj) peaks at momentum q = it 
and q < n, respectively, for points localized on the 
left (commensurate phase) and right (incommensurate 
phase) sides of the disorder line in the phase diagram. 13 

The fact that the states |0) and |0) are no longer 
degenerate implies that the low-lying excitations are 
now bound states of two spinons, i.e., spin-1 confined 
objects. 19 Indeed, a finite value of 5 introduces an ef- 
fective attractive potential between two spinons: by con- 
verting one particular singlet of the state |0) into a triplet 
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FIG. 2: Schematic ground-state phase diagram of the frus- 
trated dimerized Heisenberg model Q. For 5 = 0, the system 
is gapless for J2 < J2C = 0.2411 and spontaneously dimer- 
ized for J2 > J2C , whereas 5 > induces a gap for all J2 . On 
the so-called disorder line, 2J2 +5 = 1 (dashed), the ground 
state is given by a singlet product state on the strong bonds; 
moreover, the disorder line separates regimes of commensu- 
rate (C) and incommensurate (INC) spin correlations. The 
shaded area is the weakly confined regime of main interest in 
this paper, with £ con f > £ma g . 



and then separating the two spinons as in the state \nm), 
the energy of the system increases by an amount propor- 
tional to \m — n\5. Such an increase is due to the fact 
the singlets between the two spinons are now located at 
the weak bonds [compare Figs.[TJ(a) and (c)]i2£ As noted 
above, the spatial extension of the lowest-energy spinon- 
spinon bound stated provides a measure of the confine- 
ment length £ C onf, with £ con f ^00 for S — > 0. 

In the following, we shall be interested in the spin- 
gapped regime of weak confinement, with ^ con f S> £mag- 
Thus we shall work at 5 -C 1 and J 2 close to 1/2, Fig. [5] 
As the data analysis is much simpler for commensurate 
(instead of incommensurate) spin correlations, the regime 
to the right of the disorder line will be touched upon only 
briefly 



B. Methods 

Ground-state properties of the spin chain fl} will be 
calculated numerically for finite-size systems with the 
aid of the density-matrix renormalization group (DMRG) 
method^ using the libraries provided by the ALPS 
collaboration^ The calculations for 2N = 100... 400 
sites employ open boundary conditions, such that strong 
bonds are located at the chain ends, Fig. Q] (a). 

In addition, near the disorder line, variational wave- 
functions consisting of nearest-neighbor singlets with ad- 
ditional spinons will be utilized, similar to those in 
Ref. 0. Details can be found in Appendix fA] 
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FIG. 3: (color online) Natural logarithm of the spin-spin cor- 
relation function (Sf =11 Sf) for J 2 = 0.45 and 5 = 10" 5 , 
10" 4 , and 10" 3 . The results are from DMRG calculations 
with m = 128 states for the spin chain {TJ with 2N — 100 
sites. The curves are vertically shifted for visualization. Only 
the sites with (Sf =11 Sf) > are shown. The solid (green) 
line is an exponential fit to the data. 



III. SPIN CORRELATIONS OF UNDOPED 
SPIN CHAINS 

As a reference, we calculate the bulk spin correlation 
function (SfSj) in the ground state of the Hamiltonian 
(P) using DMRG. We restrict our attention to small S 
and the vicinity of the disorder line. 

The typical behavior of (Sf Sf) within the bulk of 
the chain for J2 = 0.45 and several values of 5 is 
shown in Fig. |on a log- linear plot. One can clearly 
see an exponential decay of the spin correlations, i.e., 
{SfSf) ~ exp (\i — j|/£ma g ), as it is expected for a spin 
gapped system. Similar results were found for J2 = 0.48. 



In Fig.|U we show the values of £ mag , obtained from an 
exponential fit to the long-distance part to the data, as a 
function of S for J2 = 0.45 and 0.48. (In the following, all 
length scales are measured in units of the lattice spacing.) 
For a fixed value of J2, £mag increases as the dimeriza- 
tion decreases and then approaches a constant value in 
the limit S — > 0. £ m ag increases as J 2 is reduced for a 
fixed value of <5. This behavior is in agreement with the 
discussions of Sec. [II] Recall that along the line (5 = 0, 
£mag diverges as J2 — > J2C from above. We note that 
on the disorder line £ m ag is minimal, but the correlations 
cannot be reasonably fitted to an exponential decay, as 
the ground state of the infinite system is given by the 
product state |0) ©■ 



IV. SPIN CHAINS WITH A NON-MAGNETIC 
IMPURITY 

Vacancy doping in quantum paramagnets can induce 
effective magnetic moments. These impurity-induced 
moments reveal themselves in a Curie-like behavior of the 
uniform susceptibility, x ^ C/T, at intermediate temper- 
atures. Remarkably, in the presence of three-dimensional 
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FIG. 4: (color online) The magnetic correlation length for 
J2 = 0.45 and J2 = 0.48 in terms of S calculated via DMRG 
(see text for the details of the calculation). The dashed lines 
are guide for the eyes. 



couplings these induced moments can order at sufficiently 
low temperatures, thus changing the spin-gapped para- 
magnetic ground state of the pure compound into a mag- 
netically long-range ordered state upon doping. 

The appearance of effective moments upon doping va- 
cancies is a property of systems with confined spinons; for 
spin 1/2 systems it is best visualized in terms of broken 
singlet bonds where one spin is replaced by a vacancy. 
The liberated spin 1/2 is confined to the vacancy at low 
energies, resulting in an effective spin- 1/2 moment i 24 ' 25 
In contrast, in host systems with deconfined spinons, no 
moments are generated by introducing vacancies. This 
theoretical picture has been supported by various nu- 
merical studies, in particular on spin chain3£ and ladder 
systemsi 27 i 28 ' 29 

In the following, we present analytical and DMRG 
results for spin chains with a single vacancy, with fo- 
cus on the spin polarization cloud associated with the 
impurity-induced magnetic moment. It has been estab- 
lished that, in the strongly confined regime, the cloud 
consists of excited triplets, with the spatial size given 
by (magj^LS^^ In contrast, in the weakly confined 
regime the cloud size is given by the extension of the 
spinon-holon bound state, i.e., £ C onf- Indeed, we shall 
find a cloud of size max(£ mag , £ C onf), with a distinct 
crossover for £ mag ~ £ C onf- The shape of the spin po- 
larization cloud allows us to distinguish the weakly from 
the strongly confined regime. 

A remark on limits is in order here: the ground state 
of the infinite system with a single vacancy is doubly de- 
generate, corresponding to the residual vacancy-induced 
moment. Applying a uniform field h u <C A, where A is 
the bulk spin gap, fully polarizes this moment and leads 
to a static staggered spin arrangement near the vacancy. 
A finite concentration of vacancies introduces another en- 
ergy scale, namely the effective (maximum) interimpurity 
interaction J e g, see also Sec. IV CI The following discus- 
sion of the polarization cloud around a vacancy assumes 
the limit h u — ► 0, keeping J c g <§; h u <C A. 
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FIG. 5: (color online) Schematic representation of (a) the 
Hamiltonian |T} with a non-magnetic impurity (open circle) 
at site i = 1 and (b) the state \p) which consists of an empty 
site at i = 1 and a spinon at site 2p separated by a singlet 
pattern. The symbols are the same as in Fig. [1] 



Variational calculations 



0.68 - 



LU 



0.64 




FIG. 6: (color online) Variational ground state energy E as a 
function of 8 (open circles) for the doped frustrated dimerized 
spin chain. The results are for points along the disorder line, 
i.e., J2 = (1 — S)/2. The solid (blue) line is a power law fit to 
the data in the region 10~ 5 < 8 < 10~ 3 (see text for details). 



Near the disorder line the low-energy states of the 
dimerized spin chain {T]) are dominated by short-range 
singlet coverings: The 5 = 1 states \mn), Fig. [He), have 
a lower energy than any S = 2 state even for a large 
separation (m — n) between the two spinons. This is the 
basis of the variational method of Ref.[la which we adopt 
here to calculate the spin correlations in the vicinity of a 
non-magnetic impurity. 

Consider the frustrated dimerized spin chain |T]) with 
strong bonds at the chain ends and introduce a vacancy 
at an odd site. For finite dimerization, 8 > 0, the system 
has no inversion symmetry at the impurity site i; mp even 
in the infinite-system limit. In this configuration, the lib- 
erated spin-1/2 stays mainly in the region i > ii mp : recall 
that the low-energy sector is formed by states made out 
of short-range singlets. The free spin can only occupy 
a site i < Zi mp if a (long) singlet is formed between the 
spins at sites i lmp — 1 and ii mp + 1, which corresponds 
to a high-energy state. To simplify the variational cal- 
culations, we shall therefore ignore the region i < ii mp , 
i.e., consider a spin chain with the vacancy placed at the 
chain end, as illustrated in Fig. EJa). Also, we focus on 
parameter values on the disorder line, J 2 — (1 — <5)/2with 
l(T 5 <5<l(r 2 . 

We restrict the variational Hilbert space to the states 
\p) , which consist of an empty site at i = 1 and a spinon 
at site 2p in a background of nearest-neighbor singlets as 
illustrated in Fig.^b). The states \p) are not orthogonal 
to each other. In fact, the overlap between two states |m) 
and \p) reads 



(m\p) = 




A long but straightforward calculation shows that the 
matrix elements of the Hamiltonian (fT]) with a vacancy 



at site i = 1 are given by 

2(m\H-E \p) = i{2 + 5J 2 

+ 3(1 - 2J 2 ) [1 + 9 + (m,p) + (m - p)6L(m,p)] 

+ 36 [2p — 1 + 8-(m,p) + (m — p)9 + (m,p)] }(m|p) 

+ (1 + 6- J 2 )(m\p - 1) + (1 - 8 - J 2 )(m\p + 1), (5) 

whereto = -(3/8){l+S)(2N), 9+(m,p) = 9(p-m-l) + 
9(m — p— 1), and 9-(m,p) = 9(p— m — 1) — 9(m — p— 1) 
with the step function 9(x) defined as in Eq. (|A8j) . 

Instead of considering the non-orthogonal basis {\p)}, 
we can define a orthonormal one^ via the Gram-Schmidt 
orthogonalization procedure, namely 

l&> = |1>, 

I0p> = -^=[b) + (l/2)b-l>], p>2. (6) 

Diagonalizing the Hamiltonian ([5]) in the basis {|</> P )} 
gives the ground-state energy E and wavefunction 
in the presence of an impurity: 

N 

|*> = £ C pI4p>- (7) 
P =i 

In Fig. we show the energy E as a function of <5 
calculated for a spin chain with 2N = 600 sites using 
the procedure described above. The data in the region 
10~ 5 < 5 < 10~ 3 , which, as we will see in the next sec- 
tion, is deep in the weakly confined regime, can be fit- 
ted by a power law do + ai<5 2 / 3 with ao = 0.6252 and 
rii = 1.6241. Such S dependence of the ground state en- 
ergy is chosen based on the analytical findings of Uhrig 
and collaborators^ who studied a similar problem via 
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FIG. 7: (color online) Variational results for the spin polariza- 
tion (Sf) around i imp = 1. J 2 = (1 — <5)/2 and 8 = 10" 3 (black 
circles), 10~ 4 (red squares), and S — 10~ 5 (green triangles). 
Top panel: linear-linear scale. Bottom panel: same data on 
a log- linear scale, showing only the sites with (Sf) > 0; the 
solid (blue) lines are fits to the data using Eq. ((9} for the 
entire data set with i > 1 and ln[(S' l z )] > —20. 

an effective one-dimensional Schrddinger equation. They 
showed that the solutions ip( x ) °f the effective equation 
can be written as ip(x) ~ Ai(x/£ + z.;), where are the 
zeros of the Airy function Ai(x)22- and £ = (3m5/2)~ 1 / 3 
is the characteristic length scale with m being the spinon 
mass, and the ground state energy E ~ S 2 ^ 3 . 

Once the coefficients C p are known, the spin polariza- 
tion (Sf) = (*|Sf |*> = Z p , m C* m C p {<f> m \S*\<f> p ) can be 
calculated with the aid of Eq. © and the expression 

(m\Sf\p) = (-iy±{e(i-2m)6(2p-i)+8(i~2p) 

x6(2m - *)] (1 - (l/2)6 m , p ) (m\p). (8) 

The results for 5 = 1CP 5 . . . I0 -3 are shown in Fig. [Jj 
on a linear-linear and log-linear scales. One can see 
that the spin polarization has a staggered profile in the 
vicinity of ii mp as is was already realized in previous 
calculations ! 19 ' 28 ' 33 By decreasing the dimerization, not 
only the position of the maximum of (Sf) moves away 
from the impurity site but also the distribution gets 
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FIG. 8: (color online) Variational £ cm (filled red triangles), 
£max (open green circles) and f i (filled black circles) as a func- 
tion of 5 for the disorder line. The solid (red) line corresponds 
to a power law fit to the data in the region 10~ 5 < S < 10 -3 
(see text for details) and the dashed line is a guide for the 
eyes. 

broader and broader. Such features indicate that the 
system approaches a deconfmed regime as the dimeriza- 
tion decreases. Indeed, at S = 0, the spin polarization 
peaks at the center of the chain. 

We now discuss possible quantitative measures of the 
size of the polarization cloud. This is not trivial, as the 
long-distance decay deviates from a pure exponential, in 
particular for small 5. The asymptotic behavior of the 
Airy function Ai(x) suggests a fit using the expression 

log(5f ) = 61 - [|t - i imp |/ei] 1+ ^ + 62 log |t - i imp \ (9) 

for i > z; mp , with parameters 61, . 3 and £1, the latter 
being the length scale corresponding to the cloud size. 
This formula, being able to interpolate between a pure 
exponential and a faster decay, depending on whether 
63 = or 63 ^ 0, nicely fits the data as is illustrated in 
Fig. [71 Alternative definitions for the polarization cloud 
size (to the right of iimp) are given by the position of 
the maximum £ max of (Sf) relative to ?i mp and by the 
center-of-mass coordinate of the cloud^ 

The variational results for £ cm , £ max , and £1 as a function 
of 5 are plotted in Fig. [5] As expected, the cloud size 
diverges as 6 — > 0, which corresponds to the MG point, 
and all measures are much larger than £ mag ~ 1. For 
instance, the data for £ cm in the region ICP 5 < 6 < 1CP 3 
can be fitted by ao + ai<5 -1 / 3 with ao = —2.3952 and 
ai = 1.7031. Again, the 6 dependence follows from the 
results derived in Ref. I20I . 

B. Numerical results 

We now use DMRG to investigate spin chains with 
a single vacancy in a parameter regime away from the 
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FIG. 9: (color online) Schematic representation of (a) the 
Hamiltonian ([T} with a non-magnetic impurity (open cir- 
cle) at site iimp and the ground-state configurations for (b) 
strongly and (c) weakly confined regimes. The symbols are 
the same as in Fig. [T] 



disorder line. We keep m = 150 states for a chain with 
2N = 400 sites and open boundary condition, and the 
vacancy is located at site i[ mp = 101, Fig. [Ha). This 
choice is due to the fact that the spin polarization cloud 
around ii, np is asymmetric (see the discussion below) . In 
the following, we only discuss the properties of lowest 
state in the Sf ot = — 1/2 sector. 

Some representative results for the spin polarization 
around the vacancy are plotted in Fig. [10] on a log-linear 
scale (only (Sf) > are shown). For the disorder line, 
the spin polarization curves are asymmetric with respect 
to iimp. On the left side of the impurity (i < ii mp ) an ex- 
ponential tail is evident while on the right side (i > ii mp ) 
the shape of the spin polarization resembles the one de- 
scribed in the last section. In fact, the results of DMRG 
and the variational approach (solid lines in Fig. [T0|) are 
in quite good agreement. The behavior for J2 = 0.48 
line is richer: Notice that for the largest value of 5, (Sf) 
is symmetric around the impurity site with exponential 
tails. As S decreases, the curves become more and more 
asymmetric, i.e., they display features qualitatively sim- 
ilar to the ones observed for the disorder line. Indeed 
for S < 10~ 2 a pure exponential decay is no longer visi- 
ble in the region i > i- lmp . Similar behavior (not shown 
here) is observed for J2 = 0.45. We have also studied 
the vacancy-induced spin polarization for J2 > 1/2: in 
this regime, the polarization displays incommensurate os- 
cillations with an envelope similar to those in Fig. [TOl 
However, the incommensurate behavior complicates nu- 
merical fits such that we refrain from a further analysis. 

We proceed as in Sec . II V Al and perform a quantitative 
analysis of the DMRG data for (Sf). Our results are 
summarized in Fig. [11] which shows the decay lengths 
£1, £2 (described below), the center-of-mass position of 
the cloud, £ cm (fit)]) , and the bulk correlation length £ m ag 
(Sec. IIIII) as function of 5 for three different values of 
Ji. We also include the variational £ cm determined in 
the previous section. One can see that the variational 
and DMRG results for £ cm along the three lines are in 
quite good agreement. This implies that the restriction 



of the Hilbert space to the subspace {\p)} is a reason- 
able approximation to describe the (weak) confinement 
physics. Recall that \p) includes only singlets made out 
of neighboring sites and the ones between the i; mp and 
the spinon site are located at the weak bonds in contrast 
to the ground state of the undoped chain whose singlet 
pairs are placed at the strong bonds. 

The spatial decay of the spin polarization cloud, (Sf), 
can be analyzed with the aid of the expression ^j. As 
(Sf) is asymmetric due to the broken inversion symmetry 
with respect to «; mp , we use Eq. ([9]) separately for the 
regions i > i lmp and i < i- lmp to extract the length scales 
£1 and £2, respectively. In both cases, we only include 
data points into the fit where (Sf) is above the noise 
floor of the numerical data, e.g., i[ mp < i < 140 in the £1 
fit of the S = 0.01 data in the top panel of Fig. [TO] The 
£1 and £2 resulting from the fit are shown in Fig. [11] as 
well. 

Let us discuss the results collected in Fig. [TT] In all 
panels, one can identify three distinct regimes: (a) For 
large S, we have £1 ~ £2 ~ £mag, i.e., both parts of 
the spin polarization cloud are determined by magnetic 
(bulk) correlations. We identify this as the strongly con- 
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FIG. 10: (color online) DMRG results for the spin polariza- 
tion (Sf) around the vacancy in the J1-J2S chain, for pa- 
rameters on the disorder line (top panel) and for J2 = 0.48 
(bottom panel). The system size is 2N = 400, and the va- 
cancy is at site i = 101. Only the sites with (Sf) > are 
shown. The solid lines in the top panel are the variational 
results calculated in Sec. IIV Al 
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FIG. 11: (color online) The characteristic length scales fi 
(open circles) and £2 (open triangles), £ cm (diamonds), and 
£mag (solid line) as function of 5 for J 2 = (1 — 5)/2, 0.48, 
0.45 from top to bottom, extracted from DMRG (see text 
for details). The variational £ cm (dashed-dotted line) is also 
shown for comparison. The dashed lines are guide to the 
eyes. Note that (i) £ma g cannot be meaningfully defined for 
J'2 = (1— S)/2, i.e., on the disorder line, and that (ii) J2 = 0.48 
and 5 > 0.04 corresponds to the incommensurate region. 



fined regime, (b) For small 5, £l » £ cm > £ 2 ~ £ m a g - 
Hence, the right part of the polarization cloud is de- 
termined by the confinement length, which is moreover 
much larger than the bulk correlation length. These are 
characteristics of the weakly confined regime, and we con- 
clude £1 = £ C onf- (c) There is a pronounced crossover 



regime between (a) and (b), e.g. 1CP 3 < S < 1CP 2 for 
J2 = 0.48. As anticipated above, £ con f diverges as 6 — > 
while £ ma g approaches a constant value - this proves the 
existence of two length scales in the regime of weak con- 
finement. 

Qualitatively, the behavior of the spin polarization 
cloud is easily understood: for strong confinement, i.e. 
large r5, we have £ C onf ~ 1 which implies that the spinon 
is close to limp. Fig. [9lb) illustrates the most proba- 
ble configuration of the ground state. Notice that the 
singlet pairs are at the strong bonds as in the ground 
state of the undoped chain. Therefore the enhancement 



of the spin correlations around 



'imp 



is mainly due to 



triplets which are excited by the coupling between the 
spin adjacent to the impurity (i.e. the spinon) and 
the environment i 27 i 28 ' 29 i 30 ' 31 The triplet physics is dom- 



inated by £ mag , hence £1 « £ 2 « £ n 



On the other 



hand, in a weakly confined regime of small 5, £ CO nf S> 1. 

, Fig. He), 



1 > 1 



imp 



The spinon travels to distant sites 
where it may also generate triplets. In contrast, config- 
urations with the spinon at sites i < i- lmp are suppressed 
as they require a long, i.e., high-energy, singlet. Hence, 
the impurity-induced spin polarization for i > ij mp prof- 
its from the travelling spinon and is determined by £ C onf 
(provided that £ C onf > (mag) whereas the polarization for 
i < « imp decays on the bulk scale £ m ag- 

To summarize, the size and shape of the spin polariza- 
tion induced by a non-magnetic impurity allow to distin- 
guish the regimes of weak and strong confinement. In the 
present ID case, both length scales £ C onf and £ m ag can in 
principle be extracted from (Sf). 



V. WEAK CONFINEMENT IN EXPERIMENTS 

In this section, we briefly sketch how weak confinement 
manifests itself in standard experiments. As before, the 
discussion is focused on ID systems, but qualitatively 
applies to D > 2 as well, see Sec. IVI1 



A. Dynamic susceptibility: Multiple bound states 

We start with the zero-temperature dynamic suscep- 
tibility x(Qi w ) °f a clean (i.e. impurity-free) gapped 
confined paramagnet. The minimum of the triplon dis- 
persion is located at some wave vector q — Q, where 
x(<f, u) — 0) displays a maximum. 

In the conventional case of strong confinement, 
x"{Q,oj) consists of a single sharp peak at the spin-gap 
energy A, corresponding to the triplon particle, and a 
continuum of multi-particle states, which typically starts 
at 3A. (In low-symmetry situations, like the J1-J2S 
chain, a two-particle continuum starting at 2A at q = 2Q 
is present.) In contrast, in a situation of weak confine- 
ment, there will be additional sharp peaks between A 
and the onset of the continuum. This can be easily un- 
derstood: while A is the energy of the lowest spin-1 two- 
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spinon bound state, weak confinement implies that higher 
two-spinon bound states are energetically only slightly 
above A. The crossover to strong confinement may be 
defined by the criterion that A becomes equal to the ener- 
getic difference between the two lowest S = 1 two-spinon 
states. 

For the J1-J2S chain considered above, the dynamic 
structure factor has been discussed e.g. in Ref. [33- The 
energy spacing of the two-spinon bound states scales as 
<5~ 2 / 3 , i.e. the number of bound states below the contin- 
uum diverges as 5 — > 0. The precise structure of x"(Q, lo) 
is determined by the energies and matrix elements of the 
two-spinon bound states, the latter giving the weights of 
the peaks in x" ■ The energies depend on the bound state 
number n as A + n 2 / 3 . Ref. |35| also provides an estimate 
of the peak weights, but the employed continuum limit 
is questionable because the weight depends on the short- 
distance behavior of the two-particle wavefunction which 
is strongly influenced by lattice effects due to the non- 
analytic nature of the confinement potential V(x) = a\x\. 
Numerical results fo r x"(q , lo) in the J1—J2—S chain can be 
found in Refs. 19.36.37.38, which clearly show additional 
sharp modes below the continuum. In particular, Fig. 4 
of Ref. [H reports energies of triplon bound states on the 
disorder line for 5 = 0.4. As both single-triplon states 
and two-triplon spin-1 bound states correspond to two- 
spinon spin-1 bound states, their results can be directly 
cast in our language: there are two (four) two-spinon 
spin-1 bound states below the two-particle continuum 
at wave vector (tt/2), but already the second is very 
close to the lower edge of the continuum, consistent with 
8 = 0.4 being located in the crossover regime between 
weak and strong confinement, Fig. 111! 
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FIG. 12: (color online) Distribution function I(x) (NMR spec- 
trum) with broadening e = 0.02 for J2 = 0.48 and several val- 
ues of 5 calculated using the DMRG results of Sec. II V Bl The 
curves are vertically shifted for better visualization. From 
bottom to top we have S = 0.0001 (black), 0.004 (red), and 
0.04 (blue). Inset: I(x) with a smaller broadening e = 0.01. 



component of the polarization, which can be defined in 
ID as^i 



(Sf W s = 7 (-l) 1 [2(Sf ) - <Sf +1 ) - (Si,)] . (11) 



For the spin chain of Sec. HVBi we calculate the distri- 
bution of staggered Knight shifts using 



27V 



i{ X ) = y j u[x-(sf) 



1=1 



stagj , 



(12) 



B. NMR line shape due to vacancies 

We now return to confined paramagnets doped with 
non-magnetic impurities, which generate a spin polariza- 
tion cloud as discussed in detail in Sec. IIVI This spin 
polarization can be made static by applying an exter- 
nal field larger than the interimpurity interaction (see 
Sec. IV CI below). The method of choice would be spin- 
resolved scanning tunneling microscopy (STM) imaging 
of the cloud which, however, is not readily available with 
atomic resolution. Alternatively, nuclear magnetic reso- 
nance provides a tool to detect the distribution of local 
fields. For instance, for doped antiferromagnetic spin 
chains^ a broadening of the NMR spectrum at very 
low temperatures is observed while, for a doped Haldane 
chain, the satellite peaks in the NMR spectrum even al- 
low the experimental determination of the magnetic cor- 
relation length^ 

The Knight shift Ki of the resonance field at the site i 
is proportional to the local electronic moment (omitting 
the factor of the electron- nucleus hyperfine coupling); the 
NMR spectrum simply corresponds to the distribution of 
the local shifts Ki. We are interested in the staggered 



where f e (x) is chosen as a Lorentzian with a small width 
e in order to have a smooth I(x) curve. 

In Fig. [12] we show examples of the NMR spectrum 
calculated using Eq. (fT2"]) and the DMRG results derived 
in Sec. II V Bl for a chain with a single impurity. The re- 
sults are for J 2 = 0.48 and S = 0.0001, 0.004, and 0.04 
which are in the regions of weak, intermediate, and strong 
confinement, respectively. As it was already observed in 
Ref. the main effect is the broadening of the spectral 
line. 

Weak confinement leads to an asymmetric spectrum, 
with a pronounced shoulder for one sign of x. This asym- 
metry is directly related to the spatial asymmetry of the 
spin polarization cloud and can be understood as follows: 
for an impurity as in Fig. [9] sites with i > i- nrip (i < z; mp ) 
have (Sf)stag < ((-Sf)stag > 0), respectively. Hence, 
the part of the cloud dominated by weak confinement, 
i > kmpi contributes to I(x < 0), where the large exten- 
sion of the cloud results in a prominent tail of I(x < 0). 
Obviously, for a finite concentration of impurities, evenly 
distributed on both sublattices, the asymmetry of the line 
shape will disappear, but the shoulder feature remains. 
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FIG. 13: (color online) Schematic representation of the Hamil- 
tonian (0 with two non-magnetic impurities (top panel) and 
two spinons between two impurities separated by a distance 
L (bottom panel). The symbols are the same as in Fig. [1] 



C. Effective interaction between vacancy-induced 
moments 

In the experimentally relevant situation of a finite con- 
centration of non-magnetic impurities, the induced mo- 
ments couple via an indirect exchange interaction me- 
diated by the gapped host. If the underlying lattice is 
bipartite, then the arising magnetism is typically non- 
frustrated, despite the strong disorder inherent to the 
problem. In D = 3 and at sufficiently low temper- 
atures, this interaction may lead to impurity-induced 
magnetic order in the paramagnetic host. Experimen- 
tally, this phenomenon has been observed in a variety 
of magnets, e.g., Zn-doped CuGe03 ) 41 i 42 i 43 Zn-doped 
SrCu 2 3 ^ and Mg-doped TlCuCl 3 4£ 

On the theoretical side, vacancy-induced magnetic or- 
dering in such compounds is commonly described by ef- 
fective models, which exclusively contain the impurity- 
induced moments of spin 1/2, which are the only low- 
energy degrees of freedom below the bulk spin gap (for 
small impurity concentration). The coupling between 
these moments is written down as a Heisenberg-type 
model, with spins Si at the random impurity locations, 
and a pairwise exchange couplings Jij of the form 



Jij oc (-l) r ' j exp(-r,j/£ mag ). 



(13) 



Here is the distance between two impurities at sites 
i and j. Note that the common assumption is that the 
decay length of the Jij is the bulk correlation length £ mag . 

In the following, we argue that Eq. (fl"3"|) does not ap- 
ply to the regime of weak confinement: there, the decay 
length of the Jij is given by £ C onf instead. The simplest 
argument is that Jij is determined by the overlap between 
the spin polarization clouds associated with the vacancies 
at sites i and j. As shown above, the cloud size crosses 
over from £ mag in the strongly confined regime to £ con f 
in the weakly confined regime - the same now applies to 
the decay length of the Jij . 

To make this explicit, let us consider the spin chain 
described by the Hamiltonian (Q]) with two vacancies, the 
first one located at an odd site and the second at an even 
site, as illustrated in Fig. [13J In this configuration, the 
two spinons stay mainly in the region of size L between 
the two impurities. As for the strongly confined regime, 
we assume that in the ground state the two spinons are 



in a singlet configuration with energy Eg =0 and define 
the effective exchange coupling as 



J' 



eff 



E S=1 



(14) 



where Eg =1 is the energy of the lowest triplet 
configuration^ 

We proceed along the lines of Refs. [2fj||35l and study the 
two-impurity problem near the disorder line via an effec- 
tive Schrodinger equation in the continuum limit, which 
captures the physics in the Hilbert space sector without 
excited triplets. The continuum limit is justified here, 
as the effective interaction turns out to be dominated 
by the long-distance tail of the wavefunction. Within 
this approach, the difference between the singlet and the 
triplet configurations is given by the boundary condi- 
tions imposed on the total wavefunction ty(x,x'). For 
the triplet configuration, the two spinons are not allowed 
to penetrate through each other while such process is 
permitted for the singlet one. Moreover, if we assume 
that the total wavefunction can be written as a prod- 
uct, i.e., ^(x, x') ~ ip(x)ip(x'), then we can reduce the 
two-particle problem to a single-particle one. Such con- 
siderations lead us to describe a spinon of mass m at 
position x by the Schrodinger equation 

subject to the following boundary conditions 



singlet: tp s (x = 0) = ip s (x = L) = 0, 



(16) 



triplet: ip^x = 0) = ^{x = L/2) = 0. (17) 



Notice that Eq. (jl6|) reflects the fact that the motion 
of the spinon is restricted to the region between the two 
impurities while the second term of Eq. (|17[) simulates the 
nonpenetration condition for the triplet configuration. 

The differential equation (|15p has an analytical solu- 
tion given in terms of a linear combination of shifted Airy 
functions Ai(x) and Bi(x), i.e.^ 

tp(x) = oAi(x/£ + x Q ) + pBi(x/£ + x ), (18) 

with the corresponding eigenvalue 

Jx 



E 



2m£ 2 



(19) 



Based on the discussions on Sec. IIV Al £ = (3to<5/2) -1 / 3 
is identified with £ con f . xq is determined from the bound- 
ary conditions (|16p and (|17p. more precisely, from the 
solutions of the following equation 



Bi(xo)Ai(Z + x Q ) - Ai(a;o)Bi(; + x ) = 0, 



(20) 



with I = L/£conf (singlet) and I = L/(2£ con f) (triplet). 
We solve Eq. pp|) numerically, calculate the correspond- 
ing effective exchange coupling as a function of L/£ con f 
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FIG. 14: (color online) Effective exchange coupling J e 
[Eq. I|14p] as a function of the renormalized impurity-impurity 
distance L/£ C onf • Inset: Natural logarithm of J cff (open green 
circles) and the corresponding linear fitting (dashed blue line) . 



according to Eq. (fl"4|) . and plot the result in Fig. (fl"4|) . 
One can see that J eff decays exponentially with the 
renormalized impurity-impurity distance L/£ con f. In- 
deed, the log-linear data in the inset of Fig. (fT4|) can 
be fitted by ln(J cff ) = o - al(L/£ conf ) with a = 3.8681 
and ai = 0.9491. 

As the calculation applies to a parameter regime of 
the spin chain where the bulk correlation length is of or- 
der unity, this shows that the effective interaction follows 
Eq. (fl"3"l) with the replacement £ mag — > £ C onf ■ 

D. Vacancy-induced order: Staggered 
magnetization 

In a situation where the effective interaction discussed 
above leads to long-range magnetic order, the stag- 
gered magnetization, being the order parameter, will be 
strongly inhomogeneous. Nevertheless, the Bragg peak 
seen in the elastic neutron scattering signal at the or- 
dering wave vector Q will be proportional to the total 
staggered magnetization of the system. For dilute impu- 
rities and neglecting fluctuation effects, this is just given 
by the number of impurities times the staggered mag- 
netization per impurity, M s t as . Note that this quantity 
can be much larger than 1/2, in contrast to the uniform 
moment associated with each impurity which equals 1/2. 

In the following, we calculate the staggered magneti- 
zation 

2N 

M stag = £(-l)%Sf>, (21) 

i=l 

for a spin chain with a single vacancy from the DMRG 
data of Sec. IIVBI The behavior of M sta g as a function of 
6 for the disorder, J 2 = 0.48 and 0.45 lines is displayed 
in Fig. [TS] One can see that |M stag | increases as we move 
into the weakly confined regime but does not diverge in 
the limit S — > 0, despite the diverging cloud size, £ C onf — ► 
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FIG. 15: (color online) The total staggered magnetization as- 
sociated with one impurity as a function of 5 for the disorder, 
J'2 = 0.48 and 0.45 lines calculated using the DMRG results 
of Sec. IIV Bl The solid lines are fitted to the data (see text 
for details). Inset: M sta g in terms of £ con f. Results for the 
disorder (solid black line), J2 = 0.48 (dashed red line), and 
J2 = 0.45 (dashed-dot blue line) lines. 

00, in this limit. Indeed, by fitting the curves with the 
expression 

-M stag (S) = c + ci cxp(-c 2 (S C3 ), 

we find that M stag (S -> 0) = 1.7275, 1.9186, and 2.2648, 
respectively for the disorder, J 2 = 0.48, and 0.45 lines. 

The finite value of M stag as £ C onf —> 00 is to be con- 
trasted with the situation of an impurity near a magnetic 
quantum critical point in 2 < D < 4, where £ mag — ► 00 
leads to M stag ^00^ Again, the difference can be un- 
derstood within the picture developed above: on the one 
hand, for £ con f — > 00 and finite £ mag , the polarization 
cloud is primarily made out of many magnetically neu- 
tral singlets, see Fig- El Indeed, in this case, it is possible 
to show analytical!}*^ that M stag approaches a constant 
as 8 — * 0. (This constant, however, depends on £ mag , see 
Fig. [151) On the other hand, for £ mag — * 00, the cloud 
consists of abundant triplets, which lead to M stag — > 00. 



VI. HIGHER DIMENSIONS 

So far, we have discussed one-dimensional spin sys- 
tems where deconfinement occurs generically in models 
with spin 1 /2 per unit cell, and confinement has to be ex- 
plicitly induced by dimerization. The situation in D > 2 
is quite different, as deconfinement is not generic. In- 
stead, it has been shown that the low-energy physics of 
many spin models can be written in terms of gauge-field 
theories. Here, the constituents are fractionalized parti- 
cles (spinons) being minimally coupled to a gauge field 
of Z2 or U(l) symmetry^^ The gauge theory may have 
both deconfining and confining phases. While the former 
display true fractionalization, the gauge-field dynamics 
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in the latter leads to confinement of objects carrying a 
gauge charge (e.g. spinons). The deconfmed phases dis- 
play topological order, characterized by the suppression 
of topological defects of the gauge field and by non-trivial 
ground state degeneracies. In contrast to ID systems, 
non-trivial quantum phase transitions between confined 
and deconfmed phases can be driven by varying micro- 
scopic model parameters. 

This discussion makes clear that the physical mecha- 
nisms for deconfinement in ID and higher dimensions are 
quite different. However, in both cases, a system can be 
in a weakly confined regime, i.e., close to deconfinement, 
with qualitatively similar observable consequences. In 
the following, we briefly sketch which of our results can 
be directly carried over to D > 2. 



A. Confined and deconfmed paramagnets 

In spin models with spin 1/2 per unit cell, quantum 
paramagnets can come in different varieties: valence- 
bond solid phases spontaneously break the lattice trans- 
lation symmetry and display gapped spin-1 excitations, 
i.e., confined spinons. In contrast, spin liquids do not 
break any symmetry of the Hamiltonian and feature 
deconfmed spin- 1/2 spinons, which can be gapless or 
gapped. In addition, spin liquids display low-lying gauge- 
field excitations: in the Z2 case, these are gapped vortices 
of the Z 2 gauge field, dubbed "visons" , whereas the U(l) 
case displays a linearly dispersing "photon" mode. These 
excitations carry zero spin, i.e., translate into singlet ex- 
citations in terms of the original spin variables. 

Weak confinement can for instance be expected (i) 
near a paramagnetic confinement-deconfincment transi- 
tion, e.g., from a VBS to a Z2 spin liquid on the trian- 
gular lattic o 46 i 47 ' 48 i 49 , and (ii) near a deconfmed critical 
point, e.g., a transition from a VBS to a Neel state on the 
square lattice. In both cases, the energy scale associated 
to confinement will be small in the VBS phase close to 
the quantum critical point. Consequently, there will be 
a large length scale £ con f , with £ con f ^> £ m ag, which char- 
acterizes the binding of spinons, which may be detected 
using nonmagnetic impurities as described above. In a 
VBS phase, £ C onf also plays the role of a characteristic 
domain wall thickness.— 

Candidates for a deconfmed VBS-Neel transition are 
the J\-Ji and J-Q square lattice spin- 1/2 Heisenberg 
models, where the antiferromagnetic order is frustrated 
either by a second-neighbor ( J2) or ring-exchange (Q) in- 
teraction. Whereas the transition in both cases may be 
weakly first order, observable finite-temperature proper- 
ties of the J-Q model have been found to be in agreement 
with predictions of deconfmed criticality. 50 

Indications for weak confinement have indeed been 
found in a numerical studyi 2 - of the frustrated J1-J2- 
J3 Heisenberg model on the square lattice, where the 
couplings Ji connecting i-th neighbor sites are all an- 
tiferromagnetic. For J3 = 0, this model displays a 



paramagnetic phase for 0.4 < J2/J1 < 0.6, which has 
been proposed to be of columnar— or plaquette^ 2 - VBS 
type. For J2 = 0, the paramagnetic phase appears to be 
larger, 0.3 < J3/J1 < 0,7 , 53 i 54 and again columnar— and 
plaquett o 55 ' 56 VBS states have been suggested. Ref. [TH 
studied the J\-J2~J , i model using exact diagonalization, 
near the line J2 + J3 — Ji/2 where the VBS order 
is likely most stable. Interestingly, these calculations 
showed clear indications of the existence of two length 
scales: while the magnetic bulk correlation length, £ mag , 
was rather short, a significantly larger scale could be 
identified in the decay of the spin polarization near a 
vacancy - consistent with our analysis, this is to be iden- 
tified with the confinement length £ C onf • For instance, for 
■Ji.2,3 = (1,0.3,0.2), £ma g m 0.7 and Cconf « 5.5, 12 These 
findings (i) suggest that weak confinement is common to 
many frustrated low-dimensional magnets, and (ii) sup- 
port that the paramagnetic phase of the J1-J2-J3 model 
is located near a deconfmed critical point. 

Among the various differences between ID and higher 
dimensions is that weak confinement in D > 2 implies 
low-lying singlet excitations, which are energetically lo- 
cated below the triplet gap; in the gauge-field language, 
these are, e.g., visons near a confinement-deconfinement 
transition. In contrast, in the J1-J2S chain there are no 
low-lying singlets, as any excitation requires the creation 
of two spinons with an energy cost of at least the triplet 
gap. 



B. Implications of linear confinement 

In the language of spinons coupled to gauge fields, a 
treatment of the weakly confined regime is difficult, as 
the gauge-field effects leading to confinement are non- 
perturbative. In the following, we simply assume that 
the gauge field induces weak linear pairwise confinement 
between spinons (as well as between spinons and holons) 
and discuss its consequences. In fact, linear confinement 
is a plausible result of VBS formation in D > 2J£ Note, 
however, that bound states can in principle exist in a 
situation of asymptotic deconfinement, see Sec. IVI Cl 

The simplest possible analysis employs a linear con- 
finement potential V^r) = ar in the continuum limit. 
Consider a confined paramagnet with a single vacancy, 
the latter introducing a spinon into the system. Assum- 
ing the spinon to be massive with quadratic dispersion 
and to be linearly bound to the vacancy, its motion is 
described by a Z?-dimensional version of Eq. (fl~5|) : 



J 

2m 



V 2 *(r) + aJrl'(r) = E$(r) 



(22) 



Here r is the radial distance from ii mp . The solution of 
the above differential equation can be written as 



*(r) = 



ip(r) exp(-iW), 
xj}{r)Y, m {e,ct>), 



D 
D 



2- 
3. 



(23) 
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where Y l m (9, 4>) is a spherical harmonic, 3 - Unfortunately, 
the resulting radial problem has no exact (analytical) so- 
lution. We then proceed, assuming that the wavefunction 
of the lowest energy vacancy-spinon bound state can be 
approximated by 

V'(r) = (r/0 /3 exp(-(r/0 /3 ), (24) 

where £ = (2ma) -1 / 3 and the coefficient (3 is variation- 
ally determined. This choice is again motivated by the 
asymptotic form of the Airy function Ai(x) [see Eq. ©]. 
It also satisfies the boundary conditions ^(0) = and 
rp(oo) = 0. By minimizing the total energy E — 
/dr**(r)iT*(r)/ / dr|*(r)| 2 , with H given by Eq. (|22) . 
it is easy to show that, for D = 2 and I = 0, we have 

J /?78 + 2- 2 - 3 /^T(2 + 3//3) 
~ 2to£ 2 2- 2 ( 1 +«/< 9 r(2 + 2/#) ' 

where T(x) is the Gamma function^ 2 - The above expres- 
sion assumes a minimum value E = 2.1967 J/(2m£ 2 ) 
for (3 = 1.1989. Similarly, for D — 3, one finds 
E = 2.5187 J/(2m£ 2 ) and (3 = 1.2682. It is worth 
mentioning that the same variational procedure ap- 
plied to the one-dimensional case gives (3 = 1.1428 
and E — 2.4429 J/(2m£ 2 ), which is a good approx- 
imation for the exact analytical solution of Eq. (f22) . 
E = 2.3381 J/(2m£ 2 ), calculated in Ref. H 

We conclude that the vacancy physics in the weakly 
confined regime in D > 2 is also dominated by the con- 
finement length: Eq. (|24[) has a maximum for r = £, 
which can be identify with £ con f. As we move deep into 
the weakly confined region (smaller a), £ increases. The 
spin polarization cloud around ii mp , whose characteristic 
length scale is £ = £ C onf, does not display a pure ex- 
ponential decay, but instead has a profile similar to the 
ID case. This also suggests that analyzing the numerical 
data in terms of an exponential decay, as done in Ref. [H, 
might be not appropriate in the weakly confined regime. 

Two-spinon bound states can be addressed similarly. 
The starting point is the D-dimcnsional version of 
Eq. (10) from Ref.H which is nothing but Eq. (22) with 
the replacements: to — > ^ = to/ 2 and E — > E — A. Here 
\x is the reduced mass, A is the spin gap (the energy of 
the lowest bound state), and r is the relative radial dis- 
tance between the two spinons. As a consequence, the 
characteristic length scale is now £ = {2^a)-^ 3 . The 
solutions of the effective Schrodinger equation, which de- 
scribes the two-spinon bound states above the spin-gap 
A, also have the form ([25)1 . We assume that Eq. ([24"]) is a 
good approximation for the radial part of the first bound- 
state wavefunction ip\{r) and variationally determine the 
parameter (3 and the energy E\ for the case 1 = 0. Once 
j3 is known, variational wavefunctions for higher bound 
states can be easily constructed just by multiplying (r) 
by a polynomial, e.g., 

Mr) = l(r/0 2 + C2i]Mr), 

Mr) = [(r/0 4 + c 3 2(r/e) 2 + c 3 i]^i(r), (25) 



respectively, for the second and third lowest-energy 
bound state above A. The coefficients dj are deter- 
mined by requiring that the states are orthogonal to each 
other. Notice that Eq. ([25]) satisfy the boundary condi- 
tion tp(0) = and they vanish for large r. The ansatz 
(125) follows from a careful analysis of the analytical solu- 
tions of the one-dimensional problem. Indeed, we applied 
the above procedure for the D = 1 case, compared the re- 
sults with the analytical solution reported in Ref. [35| and 
verified that it provides good estimates for the energy of 
the lowest two-spinon bound states. 

For D = 2 and / = 0, we find again (3 = 1.1989 and 
that the energies of the first (£a), second (E2) and third 
(E 3 ) lowest bound states are given by E t = e, + A with 
d = 2.1967, e 2 = 4.0025, and e 3 = 5.4875 J/(2/i£ 2 ). Ac- 
cording to the discussion in Sec. IV Al the strong-weak 
confinement crossover occurs when e\ ~ A. Similar con- 
siderations hold for D = 3. 

The above analysis suggests many parallels between ID 
and higher dimensions, but falls short of capturing short- 
distance features of the bound states as well as effects of 
non-trivial statistics of the spinons. 



C. Asymptotic confinement vs. short-range bound 
states 

The preceding discussion of confinement has been sim- 
plistic in that it assumed that the presence or absence of 
long-distance (asymptotic) confinement also determines 
the nature of the low-energy excitations. While this is the 
case for the one-dimensional J]— J%— 5 spin chain, which 
features a confinement potential of the type V{x) = a\x\ 
for all length scales, this is obviously not true in gen- 
eral. While the long-distance behavior of the effective 
interaction, mediated by the gauge field, decides about 
asymptotic confinement or deconfinement, it is its short- 
distance behavior which determines the existence and 
spatial size of possible bound states at low energies. 

In other words, it is conceivable that a phase in 
D > 2 spatial dimensions with asymptotic deconfine- 
ment of spinons and associated topological order displays 
short-distance spinon bound states^ In such a situation, 
the diagnostics proposed in this paper, i.e., an vacancy- 
induced magnetic moment with cloud size £ con f and mul- 
tiple sharp modes in \" below the continuum, will sig- 
nal (weak) confinement, whereas the true nature of the 
phase is deconfined. Obviously, distinguishing such a de- 
confined phase with bound states from a confined phase 
is extremely difficult. (Numerically, topological ground- 
state degeneracies may be used to detect deconfinement.) 

To our knowledge, no explicit microscopic examples for 
fractionalized spin liquids with spinon bound states ex- 
isting below the continuum are known to date^ which 
is also related to the difficulties in numerically studying 
the excitation spectra of the relevant microscopic mod- 
els. Possible candidates are states with Z2 fractionaliza- 
tion which are believed to generically display a non-zero 
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spinon pairing amplitude. 59 It should be kept in mind 
that the bound states in question, which are determined 
by short-range physics, are not necessarily part of effec- 
tive low-energy gauge-field descriptions. 

It is worth noting that a spinless impurity in the VBS 
phase near a deconfined critical point has been showr*^ 
to bind a spinon on a length scale much shorter than the 
VBS (or confinement) length. 
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APPENDIX A: VARIATIONAL APPROACH FOR 
THE DISORDER LINE 



VII. SUMMARY 

In this paper, we have argued that quantum param- 
agnets, in particular those with frustration, can display 
an interesting regime of weak confinement of spinons. 
Here, the long-distance properties are conventional, with 
elementary excitations being spin-1 triplons; however, 
the physics of spinons is visible at shorter length scales. 
Such weakly confined paramagnets are characterized by 
two potentially large length scales, namely the com- 
mon magnetic correlation length £ mag and a confinement 
length £ C onf, which measures the typical distance be- 
tween weakly bound spinons. Weak confinement implies 
Cconf £mag; the crossover to strong confinement occurs 

at £conf ~ £mag ■ 

We have illustrated the physics of the weakly confined 
regime with explicit calculations for frustrated dimerized 
spin chains, where weak confinement can be realized for 
weak explicit dimerization on top of strong spontaneous 
dimerization. In particular, we have studied the static 
response of the system to nonmagnetic impurities: the 
staggered polarization cloud induced by the impurity de- 
cays on a length scale given by max(£ con f , £ m ag)- We 
have also briefly discussed the structure of the dynamic 
susceptibility in the clean case and the effective interac- 
tion between impurity-induced moments. Qualitatively, 
our results remain valid in higher dimensions as well: in 
D > 2, weak confinement is reached near a confinement - 
deconfinement transition, e.g., between a gapped spin liq- 
uid and a valence-bond solid, or near a deconfined quan- 
tum critical point between a Neel and a VBS phase. 

However, a caveat in D > 2 is that long-distance con- 
finement (deconfinement) of spinons is not necessarily 
equivalent to the presence (absence) of short-distance 
spinon bound states: In principle, bound states forming 
conventional triplon excitations can exist even in decon- 
fined spin liquids. This tantalizing scenario is left for 
future investigations. 

Last not least, upon inclusion of mobile charge car- 
riers, weak confinement can be expected to lead to an 
almost spin-charge separated metallic state. This may 
be a conventional Fermi liquid, but with small quasipar- 
ticle weights and strong incoherent continua in spectral 
functions: The carriers will get dressed by a large cloud of 
singlet excitations, in analogy to the dressing with triplet 
excitations near a magnetic quantum critical points 



In this section we present some details of the varia- 
tional calculations discussed in Sec. IIV Al 

We start by looking at the overlap between the states 
|m) and \p). By comparing the two states one can see 
that there is a mismatch in the position of the singlets 
in the region between the sites 2m and 2p. Each mis- 
match contributes with a factor (—1/2) and the number 
of mismatches is equal to \m — n|. Thus, Eq. ((4]) follows. 

In order to calculate the matrix elements ([5]) it is con- 
venient to introduce the operator— 

/i(i,j) = i(l-4S i -S i ). (Al) 

By applying the operator h(i,j) to a normalized singlet 
combination or a product of them, we have 

h(i,j)[i,j] = 2[iJ], 
h(i,j)[k,i]\j,l] = (-l)MM, (A2) 
h(i,j)[k,i][l,j] = [k,j][i,l] + [k,i][l,j]. 

In the presence of an unpaired spin | |)j the following 
expressions hold 

h(i,j)[i,i]\l)j = (-i)|T>iM, 

h(i,j)\m,l] = (-l)M|T)i, (A3) 

fc(i,i)MIT)i = ITMU + MIT),. 

Instead of the Hamiltonian (TT]) with a nonmagnetic im- 
purity at i; mp = 1, we consider an auxiliary one H' which 
is defined as 

H' = -2H +(Ji + J 2 )(N-1). (A4) 

The above equation can be written in terms of the oper- 
ator h(i,j) (|A1[) , namely (Ji = 1) 

N N-l 

H' = (l + 5)^2h(2i-l,2i) + (l-6)^2h(2i,2i + l) 

i=2 i=l 
N-l N 

+ J 2 h ( 2i - 1, 2i + 1) + Ja h(2i, 2i + 2). 

i=2 i=l 

(A5) 

Using the properties (|A2j) and (|A3[) we can calculate 
H'\p). For instance, the contribution of the first term of 
Eq. |(A5]) reads 

N p-l 

h(2i - 1, 2i)\p) =2(N-p)\p)-\p-l)-^\i,p), 

i=2 i=2 
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The state \i,p) is quite similar to \p) apart from the fact 
that a long singlet is introduced in the vicinity of the site 
2i, i.e., 

\i,p) = [1 , 2] [3, 4] . . . [2i - 1, 2i] [2* - 2, 2* + 1] . . . 
... | t)ap...[2iV-3,2JV-2][2iV-l,2JV]. 

For the second term of Eq. (|A5|) we have 

N-l N-l 

J2h(2i,2i + l)\p)=2(p-l)\p)-\p + l)- \i,P)' 

i—1 i=p-\-l 

with 

\i,p)' = [1,2][3,4]...| T>2 P ... 

. . . [2i, 2i + l][2i - 1, 2i + 2] . . . [2N - 1, 2N}. 

Notice that \i,p) and \i,p)' are similar, but for the former 
we have i < p while for the latter, i > p. The results for 
the two remaining terms can also be expressed in terms 
of \p), \p± 1), \i,p) and \i,p)'. With the help of (jA4| it 
is possible to show that 

2(H - E )\p) = {A P 5 + 1 + 3J 2 - 28) \p) - J 2 \ Xp ) 
+(1 + S - J 2 )\ P - 1) + (1 - 8 - J 2 )|p + 1) 

p-1 N—l 
+ (l + S-2J 2 )Y / \hP) + (l-5-2J 2 ) Y, \hP)' 

i—1 i— p+1 

(A6) 



where E = -(1/4)(1 + 25 + J 2 )(2N) and the state | Xp ) 
is given by 

|X P ) - (-l)[l,2][3,4]...[2p-2,2p+2][2p-l,2p+l] 

| t>2p ■ • • [2N — 3, 2iV — 2}[2N - 1, 2N}. 

Proceeding in the same way as for Eq. one can 
show that the overlaps between \m) and the states |Xp)> 
\i,p), and \i,p)' are given by 

(m\ Xp ) = (1/2) (mb), 

(m\i,p) = — — (m\p)[8(m — p)+8(p—m — l) 

x (0(m -i) + 46(i - m - 1))] , 

(m\i,p)' = — — (m\p) [9(p — m) + 9{m — p — 1) 

x (6(i - m) + 40(m - % - 1))] , (A7) 
where the function 6(x) is defined as 

fl W = { J III ( A8 ) 
Equation © follows from Eqs. (IA"6|) - (TATj) . 



1 T. Giamarchi, Quantum Physics in One Dimension, 
Clarendon Press, Oxford, 2003. 

2 P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006). 

3 O. Heinonen, ed., Composite Fermions, World Scientific, 
Singapore (1998). 

4 S. Sachdev, in: Quantum magnetism, Lecture Notes in 
Physics Vol. 645, edited by U. Schollwock, J. Richter, D. 
J. J. Farnell, and R. A. Bishop (Springer, Berlin 2004). 

5 S. Sachdev, Nature Physics 4, 173 (2008). 

6 R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2001). 

7 A. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 

8 A. Kitaev, Ann. Phys. (N.Y.) 321, 2 (2006). 

9 G. Misguich, C. Lhuillier, B. Bernu, and C. Waldtmann, 
Phys. Rev. B 60, 1064 (1999). 

10 T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and M. 
P. A. Fisher, Phys. Rev. B 70, 144407 (2004). 

11 A. W. Sandvik, Phys. Rev. Lett. 98, 227202 (2007). 

12 D. Poilblanc, A. Lauchli, M. Mambrini, and F. Mila, Phys. 
Rev. B 73, 100403(R) (2006); D. Poilblanc, M. Mambrini, 
A. Lauchli, and F. Mila, J. Phys.: Conds. Matter 19, 
145205 (2007). 

13 R. Chitra, S. Pati, H. R. Krishnamurthy, D. Sen, and S. 



Ramasesha, Phys. Rev. B 52, 6581 (1995). 
S. Pati, R. Chitra, D. Sen, S. Ramasesha, and H. R. Kr- 
ishnamurthy, J. Phys.: Cond. Matter 9, 219 (1997). 
K. Okamoto and K. Nomura, Phys. Lett. A 169, 433 
(1992). 

I. Affleck, in Dynamical Properties of Unconventional Mag- 
netic Systems, Vol. 349 of NATO Advanced Study Institute, 
Series E: Applied Sciences, edited by A. Skjeltrop and D. 
Sherrington (Kluwer Academic, Dordrecht, 1998), pp. 123- 
131. See also |arXiv:cond-mat/9705127fr l. 
C. K. Majumdar and D. K. Ghosh, J. Math. Phys. 10, 
1388 (1969). 

B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 47, 964 
(1981). 

E. S. S0rensen, I. Affleck, D. Augier, and D. Poilblanc, 
Phys. Rev. B 58, 14 701 (1998). 

G. S. Uhrig, F. S. Schonfeld, M. Laukamp, and E. Dagotto, 
Eur. Phys. J. B 7, 67 (1999). 

Two-spinon bound states can be of S = or S = 1 char- 
acter. The lowest S = 1 two-spinon bound state is the 
conventional triplon excitation. (For the spin chain, the 
lowest 5 = bound state is the ground state itself.) In the 
limit of weak confinement, the energetic difference between 
singlet and triplet spinon bound states becomes small for 



16 



higher excited states, as the magnetic interaction between 
two distant spinons is small. 

22 S. R. White, Phys. Rev. B 48, 10345 (1993); U. Schollwock, 
Rev. Mod. Phys. 77, 259 (2005). 

23 F. Alet et al. (ALPS collaboration), J. Phys. Soc. Jpn. 
Suppl. 74, 30 (2005); F. Albuquerque et al. (ALPS col- 
laboration), J. Magn. Magn. Mater. 310, 1187 (2007); see 
http: / /alps. comp-phys.org 

24 M. Sigrist and A. Furusaki, J. Phys. Soc. Jpn. 65, 2385 
(1996). 

25 S. Sachdev and M. Vojta, Proceedings of the XIII Interna- 
tional Congress on Mathematical Physics, London, edited 
by A. Fokas et al., (International Press, Boston 2001). 

26 B. Normand and F. Mila, Phys. Rev. B 65, 104411 (2002). 

27 T. Miyazaki, M. Troyer, M. Ogata, K. Ueda, and D. Yosh- 
ioka, J. Phys. Soc. Jpn. 66, 2580 (1997). 

28 G .B. Martins, M. Laukamp, J. Riera, and E. Dagotto, 
Phys. Rev. Lett. 78, 3563 (1997); M. Laukamp, G. B. Mar- 
tins, C. Gazza, A. L. Malvezzi, E. Dagotto, P. M. Hansen, 
A. C. Lopez, and J. Riera, Phys. Rev. B 57, 10755 (1998). 

29 C. Yasuda, S. Todo, M. Matsumoto, and H. Takayama, 
Phys. Rev. B 64, 092405 (2001), and references therein. 

30 M. Vojta, C. Buragohain, and S. Sachdev, Phys. Rev. B 
61, 15152 (2000). 

31 S. Sachdev and M. Vojta, Phys. Rev. B 68, 064419 (2003). 

32 M. Abramowitz and I. A. Stegun, Handbook of Mathemat- 
ical Functions (Dover Publications, New York, 1964) . 

33 H. Fukuyama, T. Tanimoto, and M. Saito, J. Phys. Soc. 
Jpn. 65, 1182 (1996). 

34 N. Laflorencie and D. Poilblanc, Phys. Rev. Lett. 90, 
157202 (2003). 

35 T. M. R. Byrnes, M. T. Murphy, and O. P. Sushkov, Phys. 
Rev. B 60, 4057 (1999). 

36 S. Trebst, H. Monien, C. J. Hamer, W. H. Zheng, and R. 
R. P. Singh, Phys. Rev. Lett. 85, 4373 (2000). 

37 Weihong Zheng, Chris J. Hamer, Rajiv R. P. Singh, Simon 
Trebst, and Hartmut Monien, Phys. Rev. B 63, 144411 
(2001). 

38 K. P. Schmidt, C. Knetter, and G. S. Uhrig, Phys. Rev. B 
69, 104417 (2004). 

39 M. Takigawa, N. Motoyama, H. Eisaki, and S. Uchida, 
Phys. Rev. B 55, 14129 (1997). 

40 F. Tedoldi, R. Santachiara, and M. Horvatic, Phys. Rev. 
Lett. 83, 412 (1999). 

41 M. Hase, K. Uchinokura, R. J. Birgeneau, K. Hirota, and 
G. Shirane, J. Phys. Soc. Jpn. 65, 1392 (1996). 



42 M. C. Martin, M. Hase, K. Hirota, G. Shirane, Y. Sasago, 
N. Koide, and K. Uchinokura, Phys. Rev. B 56, 3173 
(1997). 

43 L. P. Regnault, J. P. Renard, G. Dhalenne, and A. 
Revcolevschi, Europhys. Lett. 32, 579 (1995). 

44 M. Azuma, Y. Fujishiro, M. Takano, M. Nohara and H. 
Takagi, Phys. Rev. B 55, R8658 (1997). 

45 A. Oosawa, T. Ono, and H. Tanaka, Phys. Rev. B 66, 
020405 (R) (2002). 

46 S. Sachdev and N. Read, Int. J. Mod. Phys. B 5, 219 
(1991). 

47 C. Xu and S. Sachdev, Phys. Rev. B 79, 064405 (2009). 

48 R. A. Jalabert and S. Sachdev, Phys. Rev. B 44, 686 
(1991); S. Sachdev and M. Vojta, J. Phys. Soc. Japan 69, 
Suppl. B, 1 (2000). 

49 R. Moessner and S. L. Sondhi, Phys. Rev. B 63, 224401 
(2001). 

50 R. G. Melko and R. K. Kaul, Phys. Rev. Lett. 100, 017203 
(2008). 

51 N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 (1989). 

52 M. E. Zhitomirsky and K. Ueda, Phys. Rev. B 54, 9007 
(1996). 

53 P. W. Leung and Ngar-wing Lam, Phys. Rev. B 53, 2213 
(1996). 

54 L. Capriotti, D. J. Scalapino, and S. R. White, Phys. Rev. 
Lett. 93, 177004 (2004). 

55 M. Mambrini, A. Lauchli, D. Poilblanc, and F. Mila, Phys. 
Rev. B 74, 144422 (2006). 

56 M. Arlego and W. Brenig, Phys. Rev. B 78, 224415 (2008). 

57 The existence of short-distance bound states in the pres- 
ence of long-distance deconfinement is well known in quan- 
tum electrodynamics: The ( — 1/r) Coulomb potential of 
two opposite charges allows for both scattering and bound 
states (e.g. the hydrogen atom). 

58 An exception appears to be the toric code model in a 
transverse field, where bound states have been seen in the 
topologically ordered phase: J. Vidal, R. Thomale, K. P. 
Schmidt, and S. Dusuel, preprint arXiv:0902.3547 

59 T. Senthil and M. P. A. Fisher, Phys. Rev. B 62, 7850 
(2000); Phys. Rev. B 63, 134521 (2001). 

60 A. Kolezhuk, S. Sachdev, R. R. Biswas, and P. Chen, Phys. 
Rev. B 74, 165114 (2006). 

61 S. Sachdev, M. Troyer, and M. Vojta, Phys. Rev. Lett. 86, 
2617 (2001). 



